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ABSTRACT 

This  report  reviews  the  result  from  airblast  measurements 
in  models  of  three  particular  types  of  underground  ammu- 
nition storage  sites.  The  use  of  simple  scaling  relation- 
ship are  shown  to  reproduce  the  blast  data  reasonably  well 
for  a wide  range  of  geometrical  parameters.  One  of  the 
models  is  geometrically  similar  to  a large  scale  test  site 
in  rock  and  direct  comparison  of  the  model  and  full  scale 
data  shows  the  blast  reducing  effects  of  wall  roughness. 
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1 . INTRODUCTION 

*K* 

This  report  is  the  fifth  in  a series  of  five  / 1—5/ 
describing  the  results  from  an  extensive  series  of  model 
tests  on  the  blast  effects  in  underground  ammunition 
storage  sites  in  case  of  accidental  explosions.  In  the 
preceeding  reports  details  were  given  of  the  test  pro- 
gramme /I/  and  measurements  were  reported  of  chamber 
pressure  /2/,  airblast  in  the  tunnel  system  of  single 
chamber  storage  sites  /3/,  and  connected  chamber  storage 
sites  / 4/ . In  the  present  report,  data  are  presented 
on  the  air  blast  in  three  particular  types  of  storage 
sites,  one  of  which  has  been  tested  in  large  scale.  In 
particular,  direct  comparison  between  the  model  and  full 
scale  data  will  show  the  importance  of  large  wall  rough- 
ness attenuation  in  the  tunnel  system  blasted  out  of 
rock. 

The  basic  scaling  relationship  employed  in  the  present 
work  are  reviewed  briefly  in  Sec.  2,  followed  by  a 
short  presentation  of  the  experimental  details  in 
Sec.  3.  The  experimental  results  are  analysed  in 
Sec.  4,  and  the  principal  results  are  finally  summa- 
rized in  Sec.  5. 

2 . THEORETICAL  CONSIDERATIONS 

2. 1 Scaling  Relationships 

As  for  the  preceeding  reports,  the  basis  for  the  pro- 
posed scaling  relationship  in  the  present  report  is 
the  familiar  Hopkinson’s  scaling  laws  / 6/ , which  were 
reviewed  in  Report  I / 1 / . To  analyse  the  present  data, 
the  scaling  relationships  proposed  in  Report  IV  /4/, 
will  be  employed.  In  the  analysis  of  the  data  in 
Sec.  4,  direct  reference  will  be  made  to  these  earlier 
results . 


* 


Numerals  inside  slashes  refer  to  appended  'Terences. 
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2 . 2 Non-Scaling  Energy  Losses 

As  v/as  discussed  in  Report  I / 1 / , the  simple  Hopkinson's 
scaling  laws  do  not  account  for  energy  dissipating 
effects  to  the  walls,  of  which  only  the  viscous  loss 
due  to  wall  roughness  is  expected  to  he  significant. 

This  effect  is  being  dealt  with  by  using  a model  origin- 
ally proposed  by  Porzel  /7/,  which  was  reviewed  in  Re- 
port I / 1 / . However,  for  the  relatively  smooth-walled 
steel  tubes  used  in  the  present  tests,  the  effects  of 
wall  roughness  is  expected  to  be  relatively  small.  This 
will  be  discussed  in  more  detail  in  Sec.  4,  and  in 
particular  the  analysis  of  the  results  from  a large 
scale  test  will  clearly  demonstrate  the  importance  of 
wall  roughness. 

3.  EXPERIMENTAL  DETAILS 

Details  of  models,  instrumentation,  and  data  reduction 
are  to  be  found  in  Report  I / 1 / , and  we  shall  there- 
fore only  summarize  the  main  features  here. 

3. 1 Models 

Fig.  3.1a  shows  the  three  models,  Model  1,  2,  and  3 
used  in  the  tests.  Model  1 represents  a situation 
with  only  one  exit.  The  length  of  the  dead  end  tunnel 
was  L/D  = 7 and  29  as  expressed  in  tunnel  diameters, 
and  the  blast  was  measured  at  various  positions  in  the 
main  passage-way  as  well  as  at  the  dead  end.  This 
model  is  similar  to  one  of  the  models  with  two  exits 
discussed  in  Report  IY. 

Model  2 in  Fig.  3.1a  is  almost  geometrically  similar 
to  a large  scale  test  site  at  Raufoss  /8/  shown  in 
Fig.  3.1b.  The  linear  scaling  factor  between  the  model 
and  prototype  was  n = 19,7. 

There  were,  however,  some  differences  between  the  model 
and  prototype:  (l)  cylindrical  storage  chamber  and 
tunnels  in  the  model  compared  to  semi-cylindrical  ones 
in  the  prototype;  (2)  a small  bend  (20°)  in  the 
prototype  tunnel , and  (4)  significant  wall 
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roughness  in  the  prototype  tunnel  blasted  out  of  rock 
compared  to  an  essentially  insignificant  wall  rough- 
ness in  the  steel  tubes  of  the  model.  Of  these  only 
(3)  is  expected  to  produce  significantly  different 
blast  propagation  in  the  tunnel  system.  This  will  be 
analysed  more  detailed  in  Sec.  4. 

Model  3 represents  a situation  where  the  branch  passage- 
ways leading  into  adjacent  chambers  are  closed  with 
blast-proof  doors  to  prevent  detonation  propagation 
from  one  storage  chamber  to  the  next.  For  this  model, 
the  blast  was  measured  in  the  tunnel  system  as  well  as 
at  the  positions  of  the  doors. 

3.2  Explosive 

The  pressed  TNT  charges  used  in  the  tests  were  suspended 
in  the  middle  of  the  detonation  chamber  and  initiated 
•with  electrical  blasting  cap  no  8.  This  has  an  equival- 
ent TNT  weight  of  1,5  - 0,5  g. 

3.3  Instrumentation  and  Data  Reduction 

To  measure  the  pressure-time  history  in  the  models, 
standard  measurement  techniques  were  used  and  the  blast 
wave  parameters  were  evaluated  employing  special  com- 
puter programmes . The  total  incertainty  in  the  peak 
pressure  was  estimated  to  be  about  10$,  whereas  the 
uncertainty  in  the  impulse  was  estimated  to  be  about 
30$.  For  further  details,  reference  Is  made  to  Re- 
port I /l/. 

4.  EXPERIMENTAL  RESULTS  AND  ANALYSIS 

Only  the  most  significant  results  will  be  presented 
in  the  following  sections  as  part  B of  this  report, 
Report  Y B,  contains  all  the  pressure- time  recordings 
and  tabulations  of  blast  wave  parameter  values. 

4. 1 Peak  Pressure  Data  Model  1 and  2 

Because  of  the  dead  end  tunnel  in  these  models 
(Fig.  3.1a),  the  pressure-time  recordings 
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showed  certain  peculiarities  which  are  exemplified  in 
Figs.  4.1a  - 4.1c  as  smooth  curve  reproductions.  The 
most  pronounced  features  here  are  the  appearences  of 
two  or  more  peaks:  One  caused  by  the  direct  blast 

wave,  and  the  others  by  the  superposition  of  the  direct 
blast  and  the  reflected  blast  from  the  dead  end  wall. 

As  indicated,  the  peak  pressures  are  defined  as  being 
the  maximum  pressures  regardless  of  whether  this  origin- 
ates from  the  direct  blast  or  the  superimposed  blast. 

Refering  to  Report  IY  /4/,  the  peak  pressure  data  in 
the  tunnel  of  a connected  chamber  storage  site  with 
two  exits  appeared  to  scale  according  to 

c? 

P = C1(Qx+/Yt)  (4.1a) 

with 

Q = explosive  charge  weight 

V.  = combined  volume  of  storage  chamber, 

branch  passage-way,  and  main  passage- 
way up  to  the  observation  point 

x = energy  distribution  factor 

*T* 

= scaling  coefficient 
C2  = scaling  exponent 

These  earlier  results  are  now  extended  to  the  dead- 
end tunnel  system  of  Model  1 and  Figs.  4. Id  and  4.1e 
show  the  definitions  of  terms  used  in  Eq.  (4.1a)  for 
L < Lp  and  L > Lp,  respectively.  Here  is  the 
length  of  dead-end  tunnel. 

Fig  4. -If  shows  the  peak  pressure  data  p+  (towards 
the  exit)  versus  Qx  /Y+.  As  may  be  seen,  the  data 

' T*  U 

appear  to  fall  on  one  universal  curve.  To  test  this 
hypothesis,  least  squares  fits  of  Eq.  (4.1a)  were 
made  to  the  three  sets  of  data.  This  produced  the 
following  results: 
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Model  1,  L^/D  =7,5: 

p+  = (11  - 2)  (Qx+/Yt)°»37  1 °'°5  (4.1b) 

Model  1,  Ld/D  =29: 

P+  = (10  ± 3)  (Qx+/Yt)°»31  1 °'°7  (4.1c) 

Model  2,  Id/D  = 4,2: 

p+  = (8,1  ± 1,0)  (Qx+/Yt)°»37  1 °'02  (4. Id) 

Thus,  within  the  stated  error  limits,  the  data 
appear  to  he  consisted  with  the  universality 
hypothesis.  An  extended  fit  combining  all  the 
data  produced  the  following  result: 

P+  = (8  ± 1)  (Qx+/Vt)0’42  1 °’35  (4. 1 e) 

4 . 2 Comparison  Model  2 and  Prototype  Data 

As  discussed  earlier  (Sec.  3.1),  Model  2 was 
geometrically  similar  to  a large  scale  test 
site  in  rock  /8/.  In  Pigs.  4.2a  - 4.2c  direct 
comparisons  are  made  between  the  model  and 
prototype  peak  pressure  data.  As  the  prototype 
data  were  obtained  from  shock  front  velocity 
measurements  /8/,  two  sets  of  data  are  given 
for  the  model:  one  set  found  directly  from  the 

pressure-time  recordings  (open  circles)  and  one 
set  also  found  from  shock  front  velocity  measure- 
ments. As  may  be  seen,  these  two  techniques 
produce  approcimately  the  same  results. 

However,  the  prototype  data  are  much  lower 
than  the  model  data.  This  is  to  be  expected  due 
to  the  much  larger  wall  roughness  in  the  proto- 
type. To  see  the  possible  effects  of  this  wall 
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roughness,  the  prototype  data  were  analysed  using 
the  model  discussed  in  Report  III: 


Po 


= C1 (Q/Vt)  2 


(4.2a) 


Y(p0)  - Y(p)  = (C3/D)(L/D)  (4.2b) 

where : 

p = fictitious  peak  pressure  without 

wall  friction 


I 


I 


I 

* 

I 


I 


p = measured  peak  pressure 

Y(pq)  = impedanse  function  of  pQ  [Eq.  (2. 3b), Report  IIJ 

Y(p)  = impedance  function  of  p [Eq.  (2 . 3b) , Report  III 

C.j  = scaling  factor 

C2  = scaling  exponent 

C,  = 2£  e = effective  wall  roughness, 

P see  Report  I 

A nonlinear  least  squares  fit  of  Eqs.(4.2a)  and 
(4.2b)  to  the  prototype  data  produced 

P0  = (3,8  ± 1,3)(Q/Vt)°»46  1 °’°5  (4.2c) 


C3  = 0,17  ± 0,03  m 


(4. 2d) 
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This  fit  reproduced  the  data  within  an  average 
scatter  of  approximately  15/°. 

It  may  be  noted  that  the  numerical  value  for 
the  effective  wall  roughness,  = 2 € e”,  is 
quite  consistent  with  direct  observations  / 9/ . 

It  is  possible  to  extract  the  experimental 
"smooth-walled"  data  from  Eq.  (4.2b),  and  the 
results  are  shown  in  Pig.  4. 2d.  This  figure  also 
includes  the  model  data. 

As  may  be  seen,  the  prototype  data  are  consitently 
lower  than  the  model  data.  It  has  not  been  possible 
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to  account  for  this  dicrepancy,  hut  there  are  two 
factors  v/hich  may  have  constributed:  (l)  the  presence 
of  relatively  large  amounts  of  dust  and  water  in  the 
prototype  installation,  and  (2)  the  use  of  quite  old 
and  loosely  packed  TUT  charges.  In  regard  to  (l), 
Swedish  model  tests  have  shown  that  the  blast  may 
be  reduced  significantly  when  the  explosive  charge 
is  in  close  contact  with  relatively  large  amount  of 
water  (comparable  to  the  explosive  weight)  / 1 0/ . 
However,  rough  estimates  indicate  that  this  effect 
alone  will  not  account  for  the  discrepancy.  As  for 
(2),  the  TNT  used  was  many  years  old  and  could  well 
have  been  less  efficient  than  ordinary  TNT.  Besides, 
the  total  charge  was  made  up  of  relatively  loosely 
packed  TNT-filled  cartridges.  Unfortunately,  no 
free-field  calibration  tests  were  performed,  and  it 
is  therefore  not  possible  to  be  more  specific  at 
this  point. 

4.3  Impulse  and  Positive  Duration  for  Model  1 and  2 

As  usual,  the  results  for  impuls  and  positive  dura- 
tion are  more  difficult  to  combine  empirically 
than  those  for  the  front  pressure.  However,  an 
attempt  will  be  made  to  analyse  the  data  according 
to  the  proposed  scaling  relationships  in  Report  IV 
/4/  for  a connected  chamber  storage  with  two  exits: 

D? 

IAk/Q  = B1(Qx+/Vt)  exp(-F5L/D)  (4.3a) 

E 

t+Ak/Q  = E-j  (Qx+/V^.)  2exp(-E3L/D)  (4.3b) 

Pigs.  4. Id  and  4.1e  show  the  definitions  of  the 
terras  used  in  these  equations. 

Following  this  scheme,  Figs.  4.3a  and  4.3b  show  the 
scaled  impulse  and  positive  duration  data  versus 
effective  loading  density.  As  a first  approximation, 
the  results  appear  to  fall  on  discrete  curves  with 
the  proposed  empirical  functional  dependence. 
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Non-linear  least  squares  fits  of  all  the  data  to 
eqs.(4.3a)  and  (4-.  3b)  produced  the  following  results: 


IAk/Q  = (50i2)(Qx+/Yt)(_0’6li0’05)exp  (-0, 05-0, 01 )L/D  (4.3c) 
t+Ak/Q  = (27-l)(Qx+/Vt)("0»97i0’08)exp  (-0, 07-0, 01 )L/D  (4.3d) 


This  appear  to  be  quite  satisfactory  considering 
the  large  uncertainties  involved  in  the  interpre- 
tation of  these  quantities. 

4 . 4 Results  - Model  3 

Fig.  4.4a  shows  typical  smoothed  reproductions 

of  the  pressure- time  recordings  in  Model  3 
(Fig.  3.1a).  The  most  pronounced  features  here 
are  the  appearences  of  rounded  peaks,  which  are 
due  to  the  complicated  diffraction  patterns  in  this 
model.  As  for  the  other  models,  the  peak  pressures 
are  defined  as  the  maximum  pressures  and  which  not 
necessarity  are  the  front  pressures. 

, Table  4.4  shows  the  peak  pressures  and  positive  dura- 
tions at  the  dead  end  walls  in  the  branch  passage-way 
(i.e.  the  positions  where  blast  doors  would  be  situated). 
These  results  are  of  course  only  valid  for  exactly  geo- 
metrically similar  full  scale  installations. 


It  may  be  noted  that  the  peak  pressures  at  MP  2 are 
only. approximately  half  the  valves  at  MP  1 . Further- 

'Z 

more,  for  the  highest  loading  density  of  50  kg/m  , 
the  pressures  are  below  100  bar  at  MP  1 for  all  the  con- 
figurations tested.  This  would  imply  that  100  bar  doors 
will  be  sufficient  for  the  prevention  of  detonation  pro- 
pagation from  one  chamber  to  adjacent  chambers. 
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Measurements  in  the  passageway  were  made  out  to  approxi- 
mately L/D  = 55,  from  the  chamber  exit  with  loading 

•z 

densities  Q/V^  =25-50  kg/nr  for  the  four  combina- 
tions of  and  l!j  , shown  in  Pig.  5.1a.  Pig.  4.4b 
shows  the  peak  pressures  versus  an  effective  loading 
density  Qx,/^+  as  defined  earlier.  The  data  were 

1 **  t 

apparently  relatively  insensitive  to  and  , and 
no  distinction  are  made  in  the  figure.  However,  the 
data  appeared  to  fall  on  two  curves  p+  and  p_  as  shown. 
Least  squares  fits  produced 

p+  = 2,7(Qx+/Vt)0’8 

and 

p_  = 5,4(Qx_/Vt)0’5 6  . 

In  the  region  of  overlap  P_<P+>  which  seems  to 
indicate  that  the  branch  passage-way  on  the  nega- 
tive side  acted  as  an  efficient  blast  trap. 

However,  due  to  the  limited  range  of  observa- 
tions it  is  hardly  possible  to  be  more  specific. 

5.  SUMMARY  AND  CONCLUSIONS 

The  methods  developed  in  this  report  for  the  pre- 
diction of  blast  wave  propagation  in  three  special 
models  may  prove  useful  in  the  planning  of  full 
scale  installations  in  rock. 

The  results  provide  a simple  and  sufficiently 
accurate  method  for  the  determination  of  the  blast 
wave  parameters  at  any  point  in  the  transport 
tunnel  for  a wide  range  of  loading  densities  and 
geometries. 

Direct  comparisons  of  the  present  results  have  been 
made  with  one  large  scale  test  with  a similar  tunnel 
system,  but  with  a much  larger  relative  wall  rough- 
ness. The  relatively  large  difference  in  the  peak 


(4.4a) 

(4.4b) 
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pressure  attenuation  for  the  model  and  large  scale 
tests,  has  been  partly  accounted  for  as  "being  due 
to  the  difference  in  wall  roughness,  "but  there  are 
other  effects  in  the  full  scale  case  which  are  not 
well  understood. 
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Table  4.4.  Blast  data  at  measuring  points  MP1  and  MP2  in 
Model  3 for  various  test  configuration,  see 
Fig.  3.1a. 


Chamber 

volume 

( cm5 ) 

Length 

, 

L ^ (cm) 

Loading 

density 

Q/Vi(kg/m5) 

Measuring 

point 

Average  peak 
pressure 

p(bar ) 

Pos.  duration 
t+(ras) 

12150 

20 

25 

MP  1 

50 

40 

12150 

20 

25 

MP  2 

26 

50 

12150 

20 

50 

MP  1 

75 

37 

12150 

20 

50 

MP  2 

38 

43 

12150 

30 

25 

MP  1 

47 

50 

12150 

30 

25 

MP  2 

26 

40 

12150 

30 

50 

MP  1 

82 

70 

12150 

30 

50 

MP  2 

45 

60 

18500 

20 

i 

25 

MP  1 

50 

50 

18500 

20 

25 

MP  2 

27 

50 
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Fig  * 3.1a. 


Model  configuration  of  the  three  sites 
discussed  in  Report  V. 

p denotes  pressure  measurements. 
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L/P  = 7,5  (-)  (at  dead  end) 
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Pig.  4.1a.  Averaged  pressure-time  recordings  for 

Model  1 (Q  = 593  g TNT)  with  L^/D  = 7,5. 
The  minus  sign  (-)  indicates  measuring 
points  in  the  dead  end  tunnel. 
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Pig.  4.1b.  Averaged  pressure-time  recordings  for 
Model  1 (Q  = 593  g TNT)  with  Ld/D  - 29. 
The  minus  sign  (-)  indicates  measuring 
points  in  the  dead  end  tunnel. 
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Pig.  4.1c.  Averaged  pressure-time  recordings  for 
Model  2 (Q  = 39  g TNT) . 
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Pig.  4. Id.  Definitions  of  torm3  used  in  the 

coaling  relationship  for  peak  pressure 
at  distances  L < where  L ^ lo  the 
length  of  tho  dead  end  tunnel. 
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Pig.  4.1©.  Definitions  og  terns  used  in  tho 
scaling  relationship  for  pealc 
pressure  at  distances  L where 

is  tho  length  of  the  dead  end 
tunnel* 


Pig.  4.  If.  Peak  pressure  versus  effective  loading 
density  for  Hodel  1.  See  Pigs.  4. la  and 
4.113  for  definition  of  x and  V,_ . 
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5?ho  solid  line  represents  a least  squares 
fit  to  all  the  data. 
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Pig.  4.2a.  Poalt  pressure  data  compared  for  Model  2 
and  prototype.  Tho  open  and  filled 
circles  arc  the  model  data  deduced  from 
pressure- time  rocoj’dings  and  shock-front 
velocity  measurements,  respectively. 


Fig.  4.2b.  Comparison  of  Ilodol  2 and  prototype  peal: 
pressure  data. 
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Pig.  4.2c.  Comparison  of  Model  2 and  prototype  peal: 
pressure  data. 
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PiC«  4. 2d.  ,,Snooth-ual!e<l,,  poalc  pressure  data  for 

prototype  compared  with  model  data.  Ehe 
straight  lino  a r opr  e 3 out  least  squares 
fits  to  the  two  sets  of  data. 
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PiG*  4* 3a.  Scaled  impulse  data  for  Model  1 and 
2 versus  effective  loading  density. 
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Fig.  4.5b.  Scaled  positive  duration  data  for  Ilodcl  1 
and  2 versus  effective  loading  density. 
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Fig.  4.4a.  Averaged  pressure-time  recordings  for 

Model  3 (Q  = 926  g TNT,  V.  = 18500  cm3 

t ' 1 

and  = 30  cm). 
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Pig.  4*4b.  Peak  pressure  versus  offectiv  loading 
density  for  Model  3. 

She  straight  linos  represent  lea3t  squares 
fits  to  th©  data  as  discussed  in  the  text. 


